Effects of Fermion Flavor on Exciton Condensation in Double Layer Systems 
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We use fermionic path integral quantum Monte Carlo to study the effects of fermion flavor on 
the physical properties of dipolar exciton condensates in double layer systems. We find that by 
including spin in the system weakens the effective inter layer interaction strength, yet this has very 
little effect on the Kosterlitz-Thouless transition temperature. We further find that, to obtain the 
correct description of screening, it is necessary to account for correlation in both the inter layer 
and intralayer interactions. We show that while the excitonic binding cannot completely surpress 
screening by additional fermion flavors, their screening effectiveness is reduced leading to a much 
higher transition temperatures than predicted with large-N analysis. 



Dipolar fermionic condensates have been a topic of 
great interest in condensed matter physics for many 
years. In recent years, double layer systems — two quan- 
tum systems separated by a thin tunnel dielectric — have 
provided a fruitful playground in which to study dipo- 
lar superfluidity both experimentally [1-4 and theoreti- 
cally [5HH]- Interest in dipolar superfluids has received 
increased due in large part to the prediction of dipo- 
lar superfluid behavior at or above room temperature 
in double layer graphene [9UTT]. This is uniquely pos- 
sible in graphene due to the symmetric linear band- 
structure and ability to sustain large carrier concentra- 
tions in two closely-spaced layers. Yet, this prediction is 
not without significant controversy. As superfluidity is 
predicted to occur in the double layer graphene system 
outside of the quantum Hall regime, additional fermion 
flavors — beyond the top or bottom layer freedom — may 
participate in the phase transition. Theoretical disagree- 
ments over the Kosterlitz-Thouless transition tempera- 
ture (Tkt) arise from differing assumptions about the 
importance of these extra flavors for screening in dipo- 
lar exciton condensates. In the works predicting a high 
transition temperature of Tkt ^ 0.1 Tp (where Tp is the 
system Fermi temperature) the fermionic degrees of free- 
dom in the system were taken to be strongly correlated in 
a gapped condensate phase and could not screen enough 
to significantly lower Tkt- Other works that predict a 
low transition temperature of Tkt ^ 10~^ Tp assume 
that screening from additional degrees of freedom add 
independenty to effectively screen out the the interlayer 
interaction [12]|13]. While experiment will be the ulti- 
mate arbiter of the value of Tkt^ many-body theoretical 
approaches will play a significant role in understanding 
the nature of the of the phase transition. 

In this Letter, we use fermionic path integral quan- 
tum Monte Carlo (PIMC) jUl [15] to elucidate the role 
of screening in exciton condensates formed in symmet- 
ric electron-hole double layer systems. We show that 
increased fermion flavor does increase the screening in 
exciton condensates formed in symmetric double layer 



system. However, the role of screening is not as dramatic 
as predicted in previous analytic work, because of strong 
excitonic pairing and other correlations. We compare 
the static correlations and dynamic response functions of 
symmetric electron-hole double layer systems with dif- 
ferent number of fermionic flavors. We use a symmetric 
model of electrons in holes in two-dimensional sheets, 
separated by a d = 0.5 nm layer of insulating Si02. 
Through analysis topological winding numbers and pair- 
correlation functions, we show that when the spin degree 
of freedom is included in our simulations, the transition 
temperature for our system drops roughly with Tp, from 
T^'j. ^ 580 K ^ 0.19 Tf to T^j. ^ 370 0.24 Tp, 
indicating that the increased number of fermion flavors 
participating in double layer systems are not completely 
screened by the by the condensate. Furthermore, we 
show that the inclusion of intralayer interactions with 
spin- 1/2 fermions, in combination with the high layer 
carrier concentrations and attractive exciton-exciton in- 
teraction, allow the system to form biexcitons. Addition- 
ally, we use dynamic density-density response function, 
collected within the PIMC framework, to demonstrate 
that when the spin degree of freedom is included, the po- 
larizability of the system actually decreases correspond- 
ing to a complicated interplay between Pauli exclusion 
forcing identical carriers to avoid one another, increasing 
screening, and that neutral excitons do not screen charge 
effectively. 

The Hamiltonian for our system includes the kinetic 
energy of the quasiparticles and Coulomb interactions. 
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All quasiparticles have equal mass of m = 0.09me, so that 
the Fermi veloctiy approximates the velocity of quasipar- 
ticle in graphene [10]. We take the dielectric constant of 
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the tunnel barrier to be = 3.9 corresponding to Si02. 
In this work we separately consider the cases of spin- 
less and spin- 1/2 quasiparticles, corresponding to differ- 
ent numbers of fermion flavors: Nf = 2 when the carriers 
only have a layer degree of freedom, and Nf = 4 when 
the carriers have both a layer and spin degree of free- 
dom. Each of the layers in our system is assumed to be 20 
nm X 20 nm with periodic boundary conditions the plane 
of the layers. We use = 40 quasiparticles in our sim- 
ulations consisting of 20 electron and 20 holes per layer. 
The number of electrons and holes corresponds to carrier 
densities of ng = n/^ = 5 x 10^^ cm~^ and places the sys- 
tem in the exciton condensate regime pHTT] . When sim- 
ulations with spin are performed, we use 10 spin-up and 
10 spin-down particles per layer. Since the Hamiltonian 
does not have any spin-dependent terms, the spin only 
enters into the calculation when considering fermion anti- 
symmetry upon particle exchange. The Coulomb interac- 
tions are handled within the pair approximation, and the 
path is discretized with a time step hAr ~ 3 Ha"^; for 
example, 256 slices are used to represent the path when 
T = 400 K. To circumvent the well-known fermion sign 
problem in our simulations, we use a ground-state fixed 
node approximation [HI [17], as used in three-dimensional 
excitonic BEC studies [IF. At each imaginary time slice, 
we require that the electron and hole coordinates lead to 
a postive Slater determinant, det |(/)(re,i — rhj)\ where 
(p = e~^/^ is the typical BCS mean- field pairing wave- 
function typically used in quantum Monte Carlo simula- 
tions of exciton condensation in symmetric electron-hole 
systems i20|. For our simulations, we have used an 
exciton pairing radius, a, of 2.3 nm but have performed 
simulations containing a = 1 nm to 10 nm and find our 
results to be insensitive to this parameter. 

In Fig. [l] (a), we show our calculated superfluid frac- 
tion as a function of the system temperature. Our cal- 
culations reveal a suppressed Tkt when the additional 
spin degree of freedom is included. Superfluid fraction 
is an excellent way to understand phase transitions in 
systems as it saturates at 1 for low temperatures and 
asymptotes to beyond the phase boundary. We denote 
the Tkt as the temperature where the superfluid frac- 
tion has dropped to 0.5 to be consistent with previous 
work ^21j . We see a clear drop in Tkt from Tkt ^ 580 K 
in the spinless case to Tkt ^ 300 K in the spin- 1/2 case, 
illustrating the effect of increased fermion flavor on the 
phase transition. In fact, this drop in transition temper- 
ature may be understood without resorting to screening 
arguments. One would expect a decrease in Tkt simply 
from reduced quantum degeneracy as more fermion fla- 
vors are added. In an ideal Bose gas, the magnitude of 
the transition temperature is determined by the condi- 
tion n\{TKTf ^ 1, where A(T) - T-^/^ ig the thermal 
deBroglie wavelength and n is the density of identical 
Bosons. As fermions pair into bosons, there must be 
at least Nf distinguishable species of bosons, decrease 
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FIG. 1. (color online) (a) Superfluid fraction for spinless 
(black circles, Nf = 2) and spin-1/2 (blue triangles, Nf = 4) 
symmetric double-layer electron-hole condensate, calculated 
from PIMC simulations of 40 electron-hole pairs. Estimates 
of Tkt are made where the superfluid fraction drops to 1/2, 
and lines are a guide to the eye. Normalized historgrams of 
topological winding distributions, which measure superfluid 
density, Eq. ([2|, for the (b) spinless and (c) spin-1/2 case. 



the density n and hence, the transition temperature Tkt 
by 1/A^/. This effect of reduced quantum degeneracy 
degeneracy can also be seen by comparing Tkt to the 
Fermi temperature (at our density, T^^ = 3086 K and 
Tp^ = 1543 K); in each case our estimated Tkt is around 
four or five times smaller than Tp. 

The calculated superfluid fraction is estimated by the 
presence of permuting paths that wind around the peri- 
odic box In two dimensions, the superfluid fraction 
is given by. 
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In Eq. ([2|, Wx and Wy are the topological wind- 
ings (which are integer multiples of the supercell di- 
mensions) corresponding to the different path configu- 
rations. Note that for excitonic condensates, we use the 
number-coupled winding, where Wx and Wy each include 
a sum over individual windings of electron and hole quasi- 
particles. 

In Fig. |2j the bottom two rows, (a) and (b) compare 
the exciton winding probabilities for the spinless (a) and 
spin-1/2 (b) for several temperatures. At T = 150 K, 
we see similar behavior for both systems, as each sys- 
tem is nearly all superfluid. Careful examination shows 
the histogram for the spinless case have extra even-odd 
structure. In an exciton, electrons pair with the holes 
and wind as a composite particle giving rise to peaks in 
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the winding histogram at even integers. Thus showing 
clear evidence the system is an exciton condensate. The 
even-odd effect is more washed out in the spin- 1/2 state, 
which we attribute to the relatively small number, ten, 
of identical fermions. Larger simulations should sharpen 
this even-odd signature of an excitonic condensate. Note 
that biexcitons would have strong peaks at multiples of 
four; we do not see evidence of a biexcitonic condensate 
in our simulations. 

Next we examine real-space static correlation among 
the quasiparticles, due to Coulomb interactions and Pauli 
exclusion. We have collected the pair-correlation func- 
tions ^(r), which are easily calculated in PIMC by bin- 
ning the equal-time pair distributions in a histogram. In 
contrast to other QMC methods, such as variational or 
diffusion Monte Carlo, there is no bias from a trail wave- 
function, thus PIMC gives essentially exact pair corre- 
lations, aside from errors associated with the fixed- node 
approximation. In Fig. |2] we show our calculated pair 
correlation functions at T = 150 K. We find the presence 
of an exchange hole in between identical fermions, which 
is larger in the Nf = A case simply because the density 
of particles of a particular flavor decreases as 1/A^/. We 
find that there is a strong, attractive, interlayer correla- 
tion with a peak at small separations that is coin- 
cident with the formation of indirectly bounds excitons 
with holes in the top layer sitting directly above the elec- 
trons in the bottom layer. This pairing of identical spin 
states in the Nf = 4 case may be an artifact of our nodal 
model, where we have built or nodal wavefunction from 
a Slater determinant of orbitals pairing like spins only. 

Interestingly, we find a small correlation hole for the 
Nf = 4 case, arising almost equally from the electron 
and hole of opposite spin: the particle in the same layer 
(with the same charge) is slightly repelled, while the par- 
ticle in the opposite layer (with the opposite charge) is 
similarly attracted. This is a clear demonstration extra 
fermion flavors correlate among each other, even when 
they do not bind as condensed excitons. This is more 
dramatically illustrated in Fig. |2] (c) and (d). In Fig. |2j 
we plot the intralayer pair correlation in 2D correspond- 
ing to gl^ and, as expected, we see a large correlation 
hole where Pauli exclusion and Coulomb repulsion forces 
other identical particles apart. However, when we include 
spin- 1/2 particles and examine the intralayer correlation 
of gl^ in Fig. [5] (c), we see a much weaker correlation 
hole. Note that we again do not see evidence of biexci- 
ton formation, which would appear as a positive correla- 
tion of like-charged, opposite-spin particles in the same 
layer 

We do not see a large drop in Tkt as predicted in 
Refs. [12] and [131 In an effort to understand this, we ex- 
plore the dynamic correlation functions, namely the po- 
larizability. In Refs. [T2l and [T3l the large- A/'j approxima- 
tion strongly resembles the random phase approximation 
(RPA), and hence the screening they find is essentially 
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FIG. 2. (color online) Radial pair correlation functions for 
(a) spinless, Nf = 2, and (b) spin-1/2, Nf — 4 cases, (a) In- 
tralayer exchange-correlation hole for the same spin species, 
g^^{Yij). (d) Intralayer correlation hole for opposite spin 
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the Lindhard function, with suppression of screening at 
very small q due a BCS-like excitonic condensate. In 
Fig. [3) we show that our calculated polarizabilities are 
considerably different. To collect the dynamic correla- 
tion functions we sample the polarization operator be- 
tween different fermion flavors 7 and 7', 
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where riq,^ = Ji e''* ''^ 
fermion flavor 7. 

In Fig. [3] we show the n^,^/(q, 0) between different 
fermion flavors normalized by the 2D density of states as 
a function of wavelength normalized by the Fermi wave 
vector at a system temperature of T = 150 K. Here we 
find the polarization response from identical particles in 
Nf = 2 and Nf = 4 track the Lindhard function very 
well in the short wavelength limit, past q > 2k f^ as one 
would expect from Fermi liquid theory. Interestingly, we 
find similar behavior for both Nf = 2 and Nf = 4 in 
response to long wavelength excitations {q 0). This 
is expected as we are clearly in the condensate regime 
in which all carriers should be paired into excitons and, 
therefore, unable to respond to perturbations. The re- 
sponse between particles in different layers or of differ- 
ent spin is negative, indicating that these correlations 
suppress screening from the independent particle, RPA 
response. The response of the exciton is quite strong, 
which the opposite-spin flavors have a weak negative re- 
sponse, Fig.jsjb). We have calculated the polarizabilities 
at temperatures both above and below T^t, and found 
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FIG. 3. (a) Interlayer and intralayer polarizabilities for Nf = 
2 versus the perturbation wavevector normalized by the Fermi 
wave vector, (b) Interlayer and intralayer polarizabilities for 
Nf — 4 normalized by the Fermi wave vector, (c) Sum of 
the interlayer and intralayer polarizabilites for Nf = 2 and 
Nf — 4 normalized by the Fermi wave vector. In each plot 
the polarizabilities are plotted against the Lindhard response 
function (grey line) in 2D and have been normalized by the 
density of states at the Fermi energy in 2D, A^o- 



very little difference between the different temperature 
traces for both the Nf = 2 and Nf = 4 cases, indicating 
that the screening behavior is dominated by pre-formed 
excitons, with exist above Tkt- 

However, Fig. |3jc) shows that the sum of inter- and 
intra-layer polarizabiHties for Nf = 4 is smaller than that 
of Nf = 2. While the large-N calculation finds enhanced 
screening with more fermion flavors Nf which grows as 
A/", but we find a weaker trend. The reason for this differ- 
ence is excitonic binding and correlation between differ- 
ent. In the large-N expansion, the starting point is RPA 
like screening. Because each flavor can screen indepen- 
dently, it has been argue that the Coulomb interaction 
because too weak to form pairs, except for weak BCS 
pairing of states near the Fermi surface when the tem- 
perature is below Tkt- In our simulations, our model 
has strong excitonic pairing. The pairing of charge into 
neutral excitons suppresses their ability to screen. This 
excitonic pairing and correlation between quasiparticles 
of different spin is not present in the large- A^/ approxi- 
mation's perturbative expansion. 

In conclusion, we have analyzed the static and dy- 
namic response functions obtained from fermionic PIMC 



quantum Monte Carlo simulations of symmetric electron- 
hole double layer systems which are spaced a distance 
of d = 0.5 nm. We find that the addition of extra 
fermion flavors to the system reduces Tkt because the 
additional flavors are not screened out by the excitons 
in the condensate thereby weakening the interlayer inter- 
action strength which drives the transition from Fermi 
liquid to exciton condensate. However, we show that 
while fermion flavors cannot be ignored the drop in Tkt 
is not as large as predicted by the large-N expansion cal- 
culations. Our analysis of polarizability shows that the 
strong excitonic pairing in our model suppresses screen- 
ing allowing a higher transition temperature. 
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